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Abstract 

Procedural noise functions are widely used in computer graphics, from off-line rendering in movie production to 
interactive video games. The ability to add complex and intricate details at low memory and authoring cost is one 
of its main attractions. This survey is motivated by the inherent importance of noise in graphics, the widespread 
use of noise in industry and the fact that many recent research developments justify the need for an up-to-date 
survey. Our goal is to provide both a valuable entry point into the field of procedural noise functions, as well 
as a comprehensive view of the field to the informed reader. In this report, we cover procedural noise functions 
in all their aspects. We outline recent advances in research on this topic, discussing and comparing recent and 
well-established methods. We first formally define procedural noise functions based on stochastic processes and 
then classify and review existing procedural noise functions. We discuss how procedural noise functions are used 
for modelling and how they are applied to surfaces. We then introduce analysis tools and apply them to evaluate 
and compare the major approaches to noise generation. We finally identify several directions for future work. 


Keywords: procedural noise function, noise, stochastic process, procedural, Perlin noise, wavelet noise, 
anisotropic noise, sparse convolution noise, Gabor noise, spot noise, surface noise, solid noise, anti-aliasing, 
filtering, stochastic modelling, procedural texture, procedural modelling, solid texture, texture synthesis, spectral 


analysis, power spectrum estimation 


ACM CCS: 13.3 [Computer Graphics]: Picture/Image Generation—I.3.7 [Computer Graphics]: Three- 
Dimensional Graphics and Realism-Colour, shading, shadowing, and texture 


1. Introduction 


Efficiently adding rich visual detail to synthetic images has 
always been one of the major challenges in computer graph- 
ics. Procedural noise is one of the most successful funda- 
mental tools used to generate such detail. Ever since the first 
image of the marble vase, presented by Perlin [Per85] (Fig- 
ure 1), Perlin noise has seen widespread use both in research 
and in industry. Noise has been used for a diverse and ex- 
tensive range of purposes in procedural texturing, including 
clouds, waves, tornadoes, rocket trails, heat ripples, inciden- 
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tal motion of animated characters and so on. It is widely used 
both in film production and video games, and is currently im- 
plemented in every major three-dimensional (3D) computer 
graphics software package, such as Autodesk 3ds Max and 
Maya, Blender, Pixar's RenderMan®, etc. 


Procedural noise has many advantages: it is typically very 
fast to evaluate, often allowing evaluation of complex and 
intricate patterns on-the-fly, and it has a very low memory 
footprint, making it an ideal candidate for compactly gener- 
ating complex visual detail. In addition, with a suitable set of 
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Figure 1: Perlin noise. (a) Perlin’s famous noise function, 
the first procedural noise function. (Figure from [Per02], 
©ACM, 2002.) (b) Perlin’s famous marble vase, one of the 
first procedural textures created using Perlin noise. (Figure 
from [Per85], CACM, 1985.) 


parameters, procedural noise can be used to easily generate a 
large number of different patterns. Finally, procedural noise 
is often randomly accessible, so that it can be evaluated inde- 
pendently at every point in constant time. This last property 
has always been a great advantage, but takes on even higher 
significance with the advent of massively parallel GPU’s and 
multicore CPU systems. 


The most recent survey on noise is in the book of Ebert 
et al. [EMP*02]. Since then there have been a multitude of re- 
cent research results in the domain, such as [CD05, BHNO7, 
GZD08, LLDD09a], as well as many others. In this survey 
we provide a unified view of both previous techniques (e.g. 
[Per85, Pea85, Lew89, PH89, vW91, Wor96, Per02]), and 
this more recent work. We also believe that recent trends in 
hardware justify the need to take a fresh look at procedural 
noise. Since 1985, compute speed has increased much faster 
than memory bandwidth. In a sense, we can now consider 
that cycles are free, in reference to the fact that most pro- 
grams in today's architectures spend a large amount of their 
time waiting for cache misses and other kinds of memory ac- 
cess. A direct consequence is that CPU-intensive algorithms 
are becoming more and more attractive; this is one of the 
main reasons that procedural methods are regaining popular- 
ity. Periodic critical re-examination of previous and recent 
methods is thus very important. 


In this survey, we attempt to provide such a critical look at 
procedural noise methods. To provide a well-founded view 
of the field, we start with both an intuitive definition of noise 
and a formal definition based on stochastic processes in Sec- 
tion 2. In this section, we also define procedural techniques 
and the different tradeoffs implied. We then provide a high- 


level review and classification of existing procedural noise 
functions in Section 3. The following two sections examine 
the important issues of modelling details with noise (Sec- 
tion 4), and that of defining noise on a surface and how to 
perform filtering (Section 5). An important part of any sur- 
vey is to analyse the strengths and weaknesses of the various 
available methods (Section 6) and to provide a comparison 
of the different tradeoffs offered by each approach, which 
we present in Section 7. We conclude in Section 8, providing 
directions we find interesting for future work. 


2. Definition of Procedural Noise Function 


In this section, we define procedural noise function, by 
defining noise both intuitively (Section 2.1) and formally 
(Section 2.2), and by defining the adjective procedural 
(Section 2.3). 


2.1. Intuitive definition of noise 


Noise is the random number generator of computer graph- 
ics. It is a random and unstructured pattern, and is useful 
wherever there is a need for a source of extensive detail that 
is nevertheless lacking in evident structure. Random patterns 
are often described in the frequency domain. Whereas in 
the spatial domain, a signal is determined by specifying the 
value for every location in space, in the frequency domain, a 
signal is determined by specifying the amplitude and phase 
for every frequency. However, for unstructured patterns, the 
phase is random and does not contribute useful information. 
Therefore, noise is often described by its power spectrum, 
which specifies the magnitude (squared) of each frequency 
and ignores the phase. This bears some similarity to how 
a chord in music is described by a set of simultaneously 
sounding notes, each with a specific frequency. A high value 
of a specific frequency in the power spectrum corresponds 
to a high contribution of the corresponding feature size in 
the spatial domain. Noise is completely characterized by its 
power spectrum, as explained in Section 2.2. Many tasks in- 
volving noise can be described as manipulations of the power 
spectrum of the noise, or spectral control. For example, mod- 
elling a noise corresponds to shaping its power spectrum, and 
filtering a noise corresponds to damping frequencies in the 
power spectrum that are too high. 


Perlin and Hoffert [PH89] gave the following definition: 
noise is an approximation to white noise band-limited to a 
single octave. White noise contains all frequencies in equal 
mixture and with random phase, so it provides the raw ma- 
terial to generate unstructured signals with any combination 
of frequencies. A band-limited power spectrum is non-zero 
only within a specific range of frequencies. So it can be used 
as a basis in the frequency domain, that is a spectral basis, 
to shape a specific desired power spectrum for modelling or 
filtering. 
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2.2. Formal definition of noise 


A more formal definition of noise will be useful in compar- 
ing and analysing different noise constructions. We will first 
recall several definitions from random processes (see, e.g. 
Papoulis and Pillai [PP02]) and Fourier analysis (see, e.g. 
Bracewell [Bra99]). 


For a discrete-valued random process y = N(x), the nth- 
order probability density function (pdf) 


fuu. Y25 <-s Yni X1, X25 6065 Xn) 


= P(N) = yi, NQ) = ys. N Qn) = y) D 


is the simultaneous probability that the noise takes on partic- 
ular values y, at n specified locations x,. The first-order pdf 
is commonly referred to as the amplitude distribution or the 
signal histogram. 


The nth-order moments are weighted averages of the cor- 
responding nth-order pdfs. The first-order moment is the 
mean 


E[N(x)] = f sioe. (2) 


The second-order moment is the expected product of the 
noise at two locations and is termed the autocorrelation or 
autocovariance (some authors define the autocovariance as 
the autocorrelation of the signal with the mean removed): 


EINEN» = [f yiyo fui, yoixi, xo) dyidy2. (3) 


A stationary random function is one whose statistics are 
invariant to a shift in the origin of the coordinate system, 
and a random function is isotropic if its statistics are also 
invariant to rotation of the coordinate system. For a stationary 
and isotropic random function the autocorrelation reduces to 
a function of a single variable, 


E[NGQG)N Q9)] = R(xi — x2)). (4) 


The autocorrelation evaluated at zero is simply the standard 
definition of variance, R(0) = E[(N(- ) — E[N( - )])*]. 
Importantly, the power spectrum of the noise is the Fourier 
transform of the autocorrelation function of the noise. 


Most existing noise functions model or approximately 
model only the first- and second-order moments, rather than 
attempting to model the full nth-order pdf. In this respect, 
noise functions are distinguished from texture synthesis algo- 
rithms (see, e.g. Wei et al. [WLKT09]) that closely reproduce 
example textures and thus necessarily reproduce their statis- 
tics. It can be seen that even the second-order pdf requires 
a lot of information to specify and manipulate. For exam- 
ple specifying an arbitrary second-order pdf for a 2D noise 
over a 4 x 4 neighbourhood, P(N(1, 1) = yy, NC, 2) = 
ym... NQ, 1) = yai... N(4, 4) = ya) involves 2566 
numbers if the values are quantized to 8 bits. However, most 
noise functions have pdfs that are jointly normal (Gaussian). 


In this case the nth-order pdf is fully and uniquely determined 
by only the first- and second-order moments, so control of 
the noise requires specifying only the desired mean and au- 
tocorrelation function or power spectrum. 


Most noises have an approximately Gaussian intensity dis- 
tribution, but for different reasons. For example, for noises 
based on frequency filtering this is because convolution im- 
plies Gaussianity [Bra99, 17], and for sparse convolution 
noises, this is because high-density shot noise implies Gaus- 
sianity [Pap71]. More generally, noise algorithms typically 
involve a weighted sum of independent pseudo-random val- 
ues. Because the pdf of a sum of random variables is the 
convolution of the pdf of the individual random variables 
[Bra99], the resulting pdf rapidly approaches Gaussian form. 


Using the preceding definitions, we take the following as 
a definition of noise: 


A noise is a stationary and normal random process. Control 
of the power spectrum is provided, either directly, or through 
the summation of a number of independent scaled instances 
of (typically band-limited) noise. 


Note that existing noise functions were not designed with 
this definition in mind; rather, this definition summarizes the 
properties of most existing noise functions. 


In summary, noise is specified through its autocorrelation 
function, or equivalently the power spectrum. Controlling a 
noise using these statistical functions is appropriate, because 
their shape is unique and specifies the character of the noise, 
whereas values of the noise itself vary randomly. The choice 
of the second-order moments is perhaps a sweet spot. These 
statistics provide a significant amount of control and have a 
well developed mathematical theory. Although modelling a 
highly structured texture by successively reproducing further 
higher-order statistics is possible, the exercise may resemble 
constructing a square wave by Fourier summation—alternate 
approaches to the goal should be considered. It is also known 
that humans have difficulty distinguishing images that differ 
only in their higher order statistics [Jul62]. Although both the 
autocorrelation and power spectrum represent the second- 
order moments, the power spectrum is the generally chosen 
representation, perhaps because it is both familiar and easily 
interpreted. 


In this survey, noise algorithms will generally be described 
for the 2D case. Generalizations to 1D, 3D and 4D are 
straightforward. The problem of defining noise at the sur- 
face of a 3D object (solid noise and surface noise) requires 
more care as discussed in Section 5. 


2.3. Definition of procedural noise 


The adjective procedural is used in computer science to dis- 
tinguish entities that are described by program code rather 


€ 2010 The Authors 


Computer Graphics Forum © 2010 The Eurographics Association and Blackwell Publishing Ltd. 


2582 A. Lagae et al./A Survey of Procedural Noise Functions 


than by data structures. Procedural techniques are code seg- 
ments or algorithms that specify some characteristic of a 
computer-generated model or effect. For example, the proce- 
dural marble texture in Figure 1 uses algorithms and mathe- 
matical functions instead of a digital photograph to define the 
colour values. We thus define a procedural noise function as 
a procedural technique for simulating and evaluating noise. 


The advantages of a procedural noise function are the 
following: 


e A procedural noise function is extremely compact, nor- 
mally requiring a few kilobytes of space compared to 
megabytes for noise images and volumes. 


e A procedural noise function is inherently continuous, 
multi-resolution, and not based on discretely sampled 
data. ^ procedural noise function can produce noise at 
any resolution desired, from an overview to extremely 
close inspection at high resolution. 


e A procedural noise function is non-periodic, filling the 
entirety of 2-D, 3-D to n-D space. In other words, it is 
unlimited in extent and can cover an arbitrary large area 
without seams and unwanted repetition. 


e A procedural noise function is parametrized, so it can 
generate a class of related noise patterns rather than being 
limited to one fixed noise pattern. The parameters control 
the power spectrum of the noise, which characterizes the 
noise pattern. 


A procedural noise function is randomly accessible. It can 
be evaluated in a constant time, regardless of the location 
of the point of evaluation, and regardless of previous eval- 
uations. This random accessibility and independent point 
evaluation make noise functions well suited to harness 
the power of multi-pipe GPUs and multicore CPUs. 


These advantages are only potential advantages. They are 
not necessarily guaranteed, but should be considered as as- 
pirations that result in the most useful procedural noise func- 
tions. 


For more information on procedural techniques, see Ebert 
et al. [EMP*02], on which this discussion is based. 


3. Overview of Procedural Noise Functions 


In this section, we give a detailed overview of procedural 
noise functions. We classify noise functions into three cat- 
egories: lattice gradient noises (Section 3.1), explicit noises 
(Section 3.2) and sparse convolution noises (Section 3.3). 
For each of these categories, we discuss a few representative 
noise functions in detail, and give an overview of related 
noise functions. We also discuss several related methods that 
do not qualify as noise functions (Section 3.4). 


3.1. Lattice gradient noises 


Lattice gradient noises generate noise by interpolating or 
convolving random values and/or gradients defined at the 
points of the integer lattice. The representative example of 
lattice gradient noises is Perlin noise. 


3.1.1. Perlin noise 


In 1985, Perlin introduced Perlin noise, his famous procedu- 
ral noise function [Per85, Per02]. 


Perlin noise determines noise at a point in space by com- 
puting a pseudo-random gradient at each of the eight near- 
est vertices on the integer cubic lattice and then doing a 
splined interpolation. The pseudo-random gradient is given 
by hashing the lattice point and using the result to choose 
a gradient. Lattice points are hashed by successive appli- 
cation of a pseudo-random permutation to the coordinates 
to de-correlate the indices into the array of pseudo-random 
unit-length gradient vectors. The set of gradients consists of 
the 12 vectors defined by the directions from the centre of 
a cube to its edges. The interpolant is a quintic polynomial, 
which ensures a continuous noise derivative. 


Since its introduction more than two decades ago, Perlin 
noise has found wide use in graphics. Perlin noise is fast 
and simple, and has continued to be the workhorse of the 
industry. 


3.1.2. Other lattice gradient noises 


Several variations, improvements, extensions and implemen- 
tations of lattice gradient noises and Perlin noise have been 
presented. 


Terminology. Ebert et al. [EMP*02] presented several in- 
stances of lattice gradient noises and a corresponding termi- 
nology. Lattice noises are defined as a noise functions based 
on the integer lattice. Value noises, gradient noises and value- 
gradient noises are defined as noise functions based on val- 
ues, gradients or both. Lattice convolution noises are defined 
as a noise functions based on convolution. We collectively 
call these noises lattice gradient noises, because the most 
well-known noise in this category, Perlin noise, is a lattice 
gradient noise. 


Other lattices. Several authors presented noise functions 
based on other lattices than the integer lattice. Wyvill and 
Novins [WN99] presented a lattice convolution noise, based 
on a more densely and evenly packed grid, inspired by sphere 
packing. Olano et al. [OHH*02] presented simplex noise, a 
Perlin-like noise based on a simplex grid. These other lattices 
lower computational complexity and eliminate undesired di- 
rectional artefacts. 


Physically-based simulations. Several authors presented 
noise functions for physically based-simulations. Perlin and 
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Neyret [PNO1] presented flow noise, a Perlin-like noise func- 
tion for generating time-varying flow textures with swirling 
and advection. Bridson et al. [BHNO7] presented curl noise, 
a Perlin-like noise function for generating time-varying in- 
compressible turbulent velocity fields. For more details about 
noise in physically-based simulations, see Bridson ef al. 
[BHNO7]. 


Better gradient noise. Kensler et al. [KKSO8] pre- 
sented better gradient noise, three mutually orthogonal 
improvements to Perlin noise. A modified hash function 
combined with a separate gradient table improves axial 
decorrelation. A different reconstruction kernel improves 
band-limitation. A projection method improves the quality 
of noise on 2D surfaces using solid noise. Note that these 
improvements apply to several lattice gradient noises. 


Hardware implementations. Several authors presented 
hardware implementations of Perlin-like noise functions. 
Hart et al. [HCK99] presented a VLSI hardware implemen- 
tation of Perlin noise. Both Hart [Har01] and Olano [Ola05] 
presented a GPU implementation of Perlin noise. Since 2003, 
noise is an integral part of the OpenGL Shading Language 
(GLSL) [Ros06]. Spjut et al. [SKB09] presented a CMOS 
hardware implementation of better gradient noise. 


3.2. Explicit noises 


Explicit noises generate noise in an explicit manner in a pre- 
process and store it. Explicit noises are not procedural noise 
functions in the strict sense, but are very relevant neverthe- 
less, which is why we cover them here. Two representative 
examples of explicit noises are wavelet noise and anisotropic 
noise. 


3.2.1. Wavelet noise 


In 2005, Cook and DeRose introduced wavelet noise [CD05]. 
Cook and DeRose observed that Perlin noise is prone to 
problems with aliasing and detail loss, because it is only 
weakly band-limited, and introduced a new noise function 
that is almost perfectly band-limited. 


In a preprocess, a tile of noise coefficients N is created. 
These coefficients represent the noise N(x) as a quadratic B- 
spline surface. This is done by creating an image R filled with 
random noise, downsampling R to create the half-size image 
R}, upsampling R? to a full size image R‘*, and subtracting 
R}? from the original R to create N. This is illustrated in 

Figure 2. The tile of noise coefficients N is thus created 
by taking R and removing the part that is representable at 
half-size. What is left is the part that is not representable 
at half-size, that is the band-limited part. The filters used in 
the downsampling and upsampling steps are obtained using 
wavelet analysis and correspond to the analysis and refine- 
ment coefficients of the uniform quadratic B-spline basis 


Figure 2: Wavelet noise generation. (a) Image R of ran- 
dom noise. (b) Half-size image R+. (c) Half-resolution image 
R‘*. (d) Noise band image N = R — RV. (Figure from 
[CD05],OACM, 2005.) 


function. The extension to more dimensions is straightfor- 
ward. 


During runtime, once the coefficients n; have been deter- 
mined, a value of N(x) for a given x can be computed using 
any evaluation method for quadratic B-splines. A small pre- 
computed volume of noise coefficients is used and space is 
tiled with that volume. 


Cook and DeRose also identified for the first time that 
sampling a 3D noise function at a 2D surface will not result 
in a band-limited texture, even if the 3D function is perfectly 
band-limited. This is discussed in more detail in Section 5. 


3.2.2. Anisotropic noise 


In 2008, Goldberg et al. introduced anisotropic noise! 
[GZD08]?. Goldberg et al. observed that existing noise 
functions only support isotropic filtering, which involves a 
tradeoff between aliasing artefacts and loss of detail, and 
presented a new noise function that supports high-quality 
anisotropic filtering. 


! When referring to the method, we will emphasize anisotropic 
noise. 


? See Lagae et al. [LZD09] for errata and clarifications. 
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The main idea of anisotropic noise is to generate noise tex- 
tures by tiling the frequency domain into oriented subbands. 
Anisotropic noise bands are not only narrowly band-limited 
in scale, but they also have a preferred orientation. The con- 
struction of anisotropic noise is based on steerable filters 
[SF95, PSOO] that partition the frequency domain. They pro- 
vide a number of properties that are crucial for noise gener- 
ation. First, each filter defines a subband that is tightly local- 
ized in scale and orientation. Secondly, the filters implement 
an invertible transform. This implies that one can exactly re- 
cover a signal from its decomposition into subbands. Finally, 
the filters are steerable in orientation. This essentially means 
that a linear interpolation of the filters can generate a filter 
with exactly the same profile, but at an intermediate orienta- 
tion. This is useful because it avoids interpolation artefacts 
when linearly blending the subbands for appropriate noise 
filtering. 


In an off-line process, noise tiles are synthesized and stored 
as discussed earlier. Each oriented subband image is packed 
into one channel of a 32-bit RGBA image, yielding four 
orientations per texture. Typically, using four or eight bands, 
that is one or two textures, leads to a good trade-off between 
storage, rendering speed, and image quality. Note that noise 
subbands are pre-computed at a single scale only. All other 
scales are generated on the fly by simply scaling the pre- 
computed textures. 


During rendering, a pixel shader computes the final noise 
value simply as a weighted sum of noise subbands at each 
pixel. By computing appropriate weighted combinations of 
the oriented subbands at each location, any desired frequency 
spectrum on the surface can be approximated. 


Goldberg et al. used anisotropic noise to obtain surface 
noise by 2D texture mapping and compensating for paramet- 
ric distortion, and for anisotropic analytic filtering. This is 
discussed in more detail in Section 5. 


3.2.3. Other explicit noises 


Two important categories of explicit noise are stochastic sub- 
division and Fourier spectral synthesis. 


Stochastic subdivision. Stochastic subdivision was intro- 
duced by Fournier et al. [FFC82], who presented the midpoint 
displacement method, a stochastic subdivision algorithm to 
generate natural irregular fractal-like objects and phenom- 
ena, such as terrain. Lewis [Lew86, Lew87] presented gen- 
eralized stochastic subdivision, a generalization of the work 
of Fournier et al. to arbitrary autocorrelation functions. 


Fourier spectral synthesis. Fourier spectral synthesis gen- 
erates a noise with a specific power spectrum by filtering 
white noise in the frequency domain (see, e.g. Bracewell 
[Bra99]). Fourier spectral synthesis was introduced in com- 
puter graphics by Anjyo [Anj88], Saupe [Sau88] and Voss 


[Vos88], who used it to generate random fractals to simulate 
natural phenomena. The mathematical texturing function of 
Gardner [Gar84] can also be seen as Fourier spectral synthe- 
sis. Fourier spectral synthesis is often used in methods for 
explicit noises, for example by van Wijk [vW91] for spot 
noise (Section 3.3.2), and by Goldberg et al. [GZD08] for 
anisotropic noise. Fourier spectral synthesis can also be use- 
ful to generate reference solutions for noise functions for 
which the expected power spectrum is known. 


3.3. Sparse convolution noises 


Sparse convolution noises generate noise as the sum of ran- 
domly positioned and weighted kernels. Three representative 
examples of are sparse convolution noise, spot noise and Ga- 
bor noise. 


3.3.1. Sparse convolution noise 


In a series of papers between 1984 and 1989, Lewis intro- 
duced sparse convolution noise [Lew84, Lew86, Lew89], 
a framework for noise functions that offers direct spectral 
control. 


The construction of sparse convolution noise is simple: an 
arbitrary kernel k is convolved with a Poisson process noise 


Y, 
No.» ff y(u, v) k(x — u, y — v) dudv. (5) 


The Poisson process consists of impulses of uncorrelated in- 
tensity a, situated at random independently chosen locations 


(Xk, yi). 


v(x, y) 2 3 aô — xi y — Ye). (6) 
k 


The Poisson process is sparse rather than being defined at 
every pixel or point in space, hence the name sparse convo- 
lution. This use of the sparse impulse noise allows some 
computational efficiency as the convolution is effectively 
splatting the amplitude-scaled kernel only at the locations 
(Xk, Yk). 


To evaluate the noise at a particular point it is necessary to 
splat only the kernels that overlap that point. This is accel- 
erated by introducing a virtual grid where the size of a grid 
cell is equal to the radius of the kernel. The evaluation then 
considers only the kernels centred in the cell containing the 
point and those in the neighbouring cells. The coordinates of 
the cell are also used to seed a random number generator for 
generating the Poisson impulses located in that cell. More 
details and improved schemes for this step are given by Wor- 
ley [Wor96] and Lagae et al. [LLDD09a]. Although Lewis 
[Lew89] describes several optimizations such as caching 
the constructed Poisson impulses under the assumption of 
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Figure 3: Sparse convolution noise. Lower left panel: win- 
dowed sample from an image of hair. Right panel: approxi- 
mate hair texture created using 2D sparse convolution using 
this kernel. 


coherent access, the sparse convolution noise is somewhat 
slower than a single octave of Perlin noise. 


Because the power spectrum of the output of a convolu- 
tion is the product of the inputs [Bra99], and the power spec- 
trum of the Poisson impulse process is constant, the power 
spectrum of the sparse convolution noise is simply a scaled 
version of that of the kernel. Direct control of the desired 
power spectrum is thus obtained simply by choosing a ker- 
nel having that spectrum. For example, a noise sharing the 
power spectrum of a sample texture can be constructed by 
using a weighted sample of the texture as a kernel (Figure 3). 
(Note that the windowing operation slightly blurs the spec- 
trum as discussed in the signal processing and filter design 
literature). A kernel with an arbitrary power spectrum can 
be constructed with the following steps: (1) generate a white 
random noise, (2) transform it to the frequency domain (be- 
cause the transform of a white noise is also white, step 1 can 
in fact be skipped), (3) filter the transformed noise with the 
desired spectral profile, (4) transform to the spatial domain 
and (5) multiply by a spatial window to produce the kernel 
(again considering standard window design issues). 


This generality in the choice of the kernel is not without 
problems however. The construction just mentioned typically 
results in kernels that do not monotonically decay away from 
the origin. Unless the density of the Poisson impulse process 
is high, valleys in the kernel are become visible as structures 
in the synthesized noise. Although this may be desirable for 
some purposes, it is objectionable in other situations and it 
violates the definition of noise as a structureless construct. 
Gabor noise (Section 3.3.3) avoids this problem while still 
providing spectral control. 


Another way of looking at the issue is in terms of phase. 
As the density of the Poisson process is increased, the phase 
structure resulting from features in the kernel is increasingly 
randomized, whereas the power spectrum of the kernel is pre- 
served. At arange of intermediate density values it is possible 
to directly synthesize noises with some textural features, as 


shown in Figure 3. However, the advent of successful tex- 
ture synthesis methods in the last decade provides a better 
approach to this problem and clarifies that noise algorithms 
are most appropriate for the random phase case. 


Although sparse convolution provided an approach to di- 
rect spectral control, it did not make any recommendation on 
which kernel to use. In the light of recent work we see that 
the implicit suggestion of allowing any kernel is in fact not 
as useful as choosing the right kernel. 


3.3.2. Spot noise 


In 1991, van Wijk introduced spot noise [VQW91], a method 
to generate stochastic textures for the visualization of scalar 
and vector fields over surfaces. Spot noise can be seen as 
an explicit form of sparse convolution noise, computed by 
scan-conversion of the spots or by Fourier spectral synthesis 
(Section 3.2.3). Although spot noise is both an explicit noise 
as well as a sparse convolution noise, it is more relevant to 
sparse convolution noises, which is why we cover it here. 


van Wijk discusses the relation between the spot and the 
texture in detail. van Wijk hinted at several important con- 
cepts which were only later introduced in the context of noise. 
For example, texture mapping on parametric surfaces, tex- 
ture synthesis over curved surfaces as an alternative to solid 
noise, and local control by variation of the spot. 


3.3.3. Gabor noise 


In 2009, Lagae et al. introduced Gabor noise [LLDD09a, 
LLDDO9b, LLD10]. Lagae et al. further developed the 
framework of sparse convolution noise by introducing the 
Gabor kernel. 


The Gabor kernel in the spatial domain, g, is the multipli- 
cation of a circular Gaussian and a 2D cosine, 


g(x,y)= Ke?) cos [27 Fo (x cos wo + y sin wo)], 
(7) 


where K and a are the magnitude and inverse width of the 
Gaussian, and Fo and wọ the frequency and orientation of the 
cosine (Figures 4(a)-(d)). The Gabor kernel in the frequency 
domain, G, is a pair of circular Gaussians, 


Gf. fy) 
K 
= exp -zh f, + Fo coso)? + (f, + Fosin ox . 
(8) 


where the Gaussians are located at the frequency with polar 
coordinates (Fo, wo), and a is the width of the Gaussians 
[Figures 4(e)-(f)]. 
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(d) (f) 


Figure 4: The Gabor kernel used in Gabor noise. (a) Gaus- 
sian. (b) Cosine. (c) Gabor kernel. (d) Gabor kernel, 3D plot. 
(e) Fourier transform of Gabor kernel. (f) Fourier transform 
of Gabor kernel, 3D plot. (Figure from [LLDD09a], GACM, 
2009.) 


Gabor noise is a sparse convolution noise with as kernel 
the Gabor kernel 


N(x. y) = Y 1 wig(K;. ai, Foi, 995:X — Xi y — yi), (9) 


i 


where {w;} are the random weights, g is the Gabor kernel, 
and {(x;, y;)} are the random positions. Depending on how 
the parameters {K;}, (aj), (Fo;) and {@o;} vary for dif- 
ferent kernels, different kinds of Gabor noise are obtained. 
When the parameters are fixed, the power spectrum of the 
noise is that of the Gabor kernel, and an anisotropic band- 
limited noise is obtained, where wo, F and a control the 
orientation, frequency and bandwidth of the noise. When the 
parameters are varied, the power spectrum of the noise is 
that of the Gabor kernel integrated over the parameters. For 
example, when {@,;} is uniformly distributed over [0, 27r), 
an isotropic band-limited noise is obtained, where Fo and a 
control the frequency and bandwidth of the noise. Lagae et 
al. use graphical user interface widgets to specify the power 
spectrum of the noise by specifying how the parameters vary 
(Figure 5). 


Lagae et al. used Gabor noise for setup-free surface noise 
and analytic anisotropic filtering of noise. This is discussed 
in more detail in Section 5. 


3.3.4. Other sparse convolution noises 


Several extensions and implementations of sparse convolu- 
tion noises have been presented. 


Figure 5: Gabor noise. Several Gabor noise patterns. The 
top row shows the Gabor noise patterns, the bottom row 
shows the corresponding widgets. (Figure from [LLDD09a], 
OACM, 2009.) 


Shaped point processes. Lewis [Lew86] presented shaped 
point processes, one of the works which would eventually 
lead to sparse convolution noise. This work hinted at sev- 
eral important concepts that would only later be fully de- 
veloped. For example, a bandpass kernel resembling the 
Gabor kernel (as in Gabor noise), filtering of noise (Sec- 
tion 5), and spatially varying noise (as in Lagae et al. 
[LLD10]). 


GPU implementations. Sparse convolution noise can be 
implemented on the GPU using splatting (point rendering, 
scan conversion) or procedurally (using a shader). Frisvad 
and Wyvill [FW07] presented a GPU implementation based 
on point rendering of sparse convolution noise with a cubic 
kernel. Lagae et al [LLDD09a] presented a procedural GPU 
implementation of Gabor noise. 


NPR Gabor Noise. Benard et al. [BLV*10] recently pre- 
sented NPR Gabor Noise, a variant of Gabor noise for co- 
herent stylization in non-photorealistic rendering (NPR). To 
preserve both the 2D aspect of the noise and the 3D mo- 
tion of the scene, the noise parameters are defined in 2D 
screen space and the point distribution is defined in 3D object 
space. 


3.4. Related methods 


Several methods have been presented that are not procedural 
noise functions but are nevertheless highly related to proce- 
dural noise functions. 


3.4.1. Fractals 


Fractals [Man82] triggered rapid and great progresses 
of procedural noise functions. For more informa- 
tion on how to use procedural noise functions in 
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practice, see, for example the work of Musgrave 
[MKM89] and Szeliski and Terzopoulos [ST89] on fractal 
terrains. 


3.4.2. Texture basis functions 


Texture basis functions are defined as functions to generate 
patterns that can be used as a basis for generating textures. 
The most well-known texture basis function is probably the 
one of Worley. Worley [Wor96] presented a cellular texture 
basis function, a texture basis function based on distances to 
feature points randomly scattered in space, which is good for 
creating textures such as flagstone-like tiled areas, organic 
crusty skin, crumpled paper, ice, rock, mountain ranges and 
craters. The implementation of Worley’s cellular texture basis 
function is very similar to that of sparse convolution noise. 
There are several other methods that could also qualify as a 
texture basis functions, for example the method presented by 
Tzeng and Wei [TW08] for parallel white noise generation 
on the GPU. 


3.4.3. Object distribution functions 


We define object distribution functions as functions to gen- 
erate patterns that consist of objects distributed over a back- 
ground. Lefebvre and Neyret [LNO3] presented pattern- 
based procedural textures, a method to generate procedural 
textures composed of randomly distributed objects on the 
GPU. Lagae and Dutré [LD05] presented a procedural ob- 
ject function, a texture basis function for objects distributed 
according to a Poisson disk distribution, which is good for 
creating textures such as polka dots. The tile-based methods 
used in this method [Lag09] are also useful in the context 
of procedural noise functions, for example for noise tiles 
[CD05, YL08]. 


4. Modelling with Procedural Noise Functions 


Creating visually rich and interesting content from noise is 
not an easy task, essentially because the random nature of 
noise makes it difficult to control and predict the result. In 
addition, noise is often only the first component in a long 
chain of operations to achieve the end result. Most systems 
for modelling with noise are based on the concept of block 
shaders [AW90], in which a texture is described as a network 
of modules. 


We mainly focus on the design of the noise patterns them- 
selves, and refer the reader to the book Texturing & Model- 
ing: A Procedural Approach [EMP*02] for a more detailed 
account of the most useful approaches to generate terrains, 
shapes and textures from noise and procedures. 


In this section, we describe spectral control of noise (Sec- 
tion 4.1), direct editing of noise values (Section 4.2), and 
noise by example (Section 4.3). 


4.1. Spectral control of noise 


As explained in Section 2.1, noise patterns are best described 
in terms of frequency content, through their power spec- 
trum. Controlling a noise pattern through its spectrum re- 
quires some training, but is convenient once the link be- 
tween the spectrum and the visual aspect of the noise is 
understood. 


We describe next the most common approaches for spec- 
tral noise control. The first approach consists of summing 
weighted layers of band-limited noise. The second approach 
discusses the specific case of sparse convolution noises, 
which are controlled through the choice of kernel. 


We would like to note that in previous work the term band- 
limited is often used where the term band-pass would be 
more appropriate. Note that a band-limited power spectrum 
is zero beyond a specific frequency, whereas a band-pass 
power spectrum is zero outside of a frequency interval (see, 
e.g. Bracewell [Bra99] or Papoulis and Pillai [PP02]). In the 
preceding sections we use band-limited for consistency with 
previous work, but in the following sections we will use the 
appropriate term. 


4.1.1. Weighted sum of band-pass noises 


Most procedural noise functions directly produce band-pass 
noises. Each noise band corresponds to an elementary ran- 
dom pattern, with a frequency content limited to a specific 
range. Note that bands at different frequencies are easily 
obtained by scaling an initial band-pass noise. 


The very reason for which procedural noise functions are 
designed to produce band-pass noises is to let complex pat- 
terns be defined by adding several bands of noise. Each band 
is multiplied by a weight controlling its contribution to the 
final result. This idea was introduced by Perlin [Per85]. The 
final pattern is obtained as 


5 w; N(2'x), (10) 


where N is a band-pass noise function and w; is the weight 
of band 7. Successive noise layers have a principal frequency 
related by a factor of two, which is why they are often called 
octaves. Perlin initially described a noise with 1/ f spectral 
content with weights computed as 1 /2'. However, in a typical 
noise modelling tool the weights are directly exposed to the 
user, as shown in Figure 6. The spectrum of the resulting 
noise pattern is obtained as the weighted sum of the band 
spectra. 
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Figure 6: Spectral control with wavelet noise. Three- 
dimensional noise patterns with 12 bands with a Gaussian 
distribution, and 8 bands with a white distribution. The blue 
bars are the band weights. These weights can be exposed to 
the user. (Figure from [CD05], CACM, 2005.) 


Noise bands that are band-pass have little overlap in 
the frequency domain and can be seen as a spectral ba- 
sis, defining a space of noise patterns. Note however that in 
the basis analogy, only positive weights are effective. More 
specifically, it is not possible to cancel energy from a fre- 
quency band because the noise has random phase. Thus, a 
resulting noise spectrum that contains no energy at some fre- 
quencies can only be produced if the primitive noise function 
is band-pass rather than merely band-limited. 


Note that the weights do not have to remain constant in 
space: By using different weights in different locations one 
can generate patterns smoothly transitioning between differ- 
ent aspects. This fact is exploited by Goldberg et al.[GZD08] 
to cancel mapping distortions and dynamically adapt the 
noise to viewing conditions (Sections 5.1 and 5.2). 


4.1.2. Sparse convolution noises 


Sparse convolution noises are controlled through the choice 
of the kernel, because the noise has the spectrum of the kernel 
(see Section 3.3.1). This choice may vary spatially to obtain 
different appearances in different areas [VW91 ]. 


Sparse convolution noises can produce band-pass noises 
with the appropriate kernel. They are thus compatible with 
the approach of summing noise bands. However, they can 
also be used to directly generate a noise with a specific spec- 
trum, provided that a kernel having this particular spectrum 
is available. This is the case for sparse convolution noise, 
which directly produces a noise with the desired spectrum 
as illustrated in Figure 7(b). Gabor noise [LLDD09a] uses 
a kernel which can itself be controlled through a number 
of parameters. These are described through widgets directly 
manipulated by the user. These parameters give direct con- 
trol over the spectrum generated by the noise, without hav- 
ing to change the kernel. Noise can evolve from anisotropic 


band-pass patterns to more elaborate patterns, as illustrated 
in Figure 5. 


4.2. Editing noise values 


In addition to spectral control, other techniques investigate 
how to control the noise values in the spatial domain. This 
is challenging to achieve without destroying the properties 
of the noise. Lewis [Lew87] generated a noise in a mul- 
tiresolution coarse-to-fine scheme, with the fine scale values 
condition on previously specified values at the coarse scale. 
However, some of the coarse scale values can be directly 
specified by the user as shown in [Lew87, Figure 5]. Yoon 
et al.[L YLC04, YLO08] let the user directly specify a few val- 
ues of a noise field. New random numbers are generated 
ensuring that the user constraints are satisfied and that the 
noise keeps its properties (value distribution, non-periodicity, 
band-pass). 


4.3. Noise by example 


To avoid manual noise design, several authors have focused 
on finding parameters from an image. This is, however, an 
extremely challenging problem. To the best of our knowl- 
edge no satisfactory solution exists for the general case of 
procedural noises as defined in Section 2. It is important to 
note that neighbourhood-based texture synthesis approaches 
as surveyed by Wei et al. [WLKTO09] do not fall in this cat- 
egory. Consequently, this is also an exciting area of further 
research. 


Several interesting solutions exist for subclasses of tex- 
tures. Ghazanfarpour and Dischler [GD95] express the noise 
as a sum of sine waves, similarly to Gardner [Gar85]. They 
select the set of sine waves from an example image, by thresh- 
olding the magnitude of its Fourier transform. This 2D func- 
tion can then be extended to define a solid (3D) noise. The 
method is further refined in subsequent work [GD96], to 
support different aspects along different directions of a solid 
noise. In [DG97], the authors focus on geometric textures. 
They analyse 1D noise profiles and automatically generate 
procedures for them. These are then extended to 2D and 3D. 
The analysis step identifies main frequencies but also per- 
forms histogram matching between the example and the gen- 
erated noise. These spectral approaches work best when the 
textures contain strong periodicities, with clearly identified 
features in the power spectrum. Lagae et al.[LVLD10] auto- 
matically compute weights of a sum of band-pass isotropic 
noise octaves, to produce an image closely resembling an 
example. The method produces results close to early by- 
example texture synthesis approaches [HB95], with the cru- 
cial difference that the result is a procedure and can be 
efficiently point-sampled. Nevertheless, this approach can- 
not faithfully reproduce structured or anisotropic patterns. 
In contrast, Galerne et al. [GGM09] randomize the phase 
spectrum of a given texture to obtain a homogeneous and 


© 2010 The Authors 


Computer Graphics Forum © 2010 The Eurographics Association and Blackwell Publishing Ltd. 


2589 


A. Lagae et al./A Survey of Procedural Noise Functions 


(6002 ‘WOV® 
'|v60Qq'T1] wo pasvq an814) ‘ajouijsa umajoads samod pagpa &q]pipuy (ç) ‘ajpwysa umaj22ds samog (p) ‘Wpssopolsad (£) uoumquisip apnyyduy (z) ‘asion (T) ‘asiou 
40q09 (f) ‘asiou ju21pv48 sayjag (2) ‘asiou 21d04jostuy (p) ‘asiou J2]2ADM (2) ‘asiou uoimjoauoo 2s4vdg (q) ‘asiou utaq (D) 'suouounf asiou ppanpa204d fo sis&qpuy. :L mZ 


"esiou 10qer) (J) "esiou “BIS jog (9) 'esiou ‘donostuy (p) "esiou jo[9AeAA (2) 'esiou *oAuoo ‘dg (q) ‘astou Ul] Jog (€) 


ca (S) 
(v) 
© 
10 
DOSE TT : a meu UTOR IP 
i (D 


£ 


fa n us - Be Bt mS Sars ee 


ioe 


> nevy] PAETA 
BAY Groa Pet AR rv wi ra yt PS IE GE 


us 


© 2010 The Authors 


Computer Graphics Forum © 2010 The Eurographics Association and Blackwell Publishing Ltd. 


2590 A. Lagae et al./A Survey of Procedural Noise Functions 


featureless noise having the same power spectrum. Recently, 
Gilet et al. [GDS10] did a first attempt at extending the 
work of Dischler and Ghazanfarpour [DG97] and Lagae et al. 
[LVLD10] to anisotropic textures using Gabor noise. Gilet 
and Dischler [GD10] also proposed a 3D extension of the 
work of Bourque and Dudek [BD04] using a mixture of 
Gabor noise and cellular noise. 


Other approaches focus on setting the parameters of exist- 
ing procedural shaders from an example image: The goal 
is to make the shader produce an image resembling the 
example as closely as possible. Although these techniques 
do not primarily target noise patterns, they could be use- 
ful to automatically select parameters of a noise function. 
Bourque and Dudek [BD04] aim at a more generic approach, 
searching for closest matches in a database of images gener- 
ated by sampling the parameter space of many shaders. Qin 
et al.[QY02] similarly optimize shader parameters using a 
genetic algorithm. 


5. Procedural Noise Functions on Surfaces 


Noise in computer graphics is especially useful to add vi- 
sual details in renderings, through texturing. A texture ob- 
tained from a noise pattern inherits all its advantages: Non- 
periodicity, low memory cost, resolution and efficient random 
access. 


In this section, we discuss how noise patterns are typically 
mapped on surfaces (Section 5.1) as well as the challenges 
this creates for anti-aliased rendering (Section 5.2). 


Using noise to texture surfaces introduces two different, 
albeit closely related, challenges. A first difficulty is to find 
an appropriate mapping of the noise to the surface, while pre- 
serving the properties of the noise (frequency content, conti- 
nuity). A second difficulty is to adapt the noise to the view- 
ing conditions. Indeed, a noise with high frequency quickly 
produces disturbing aliasing artefacts when mapped onto a 
surface seen at an angle or in the distance. 


5.1. Noise on surfaces 


There are three methods for obtaining noise on a sur- 
face: mapping a 2D noise onto the surface using a planar 
parametrization, sampling a solid noise, or defining a noise 
directly on the surface. We refer to this latter case as surface 
noise. Note that a surface noise should retain the properties 
it exhibits in 2D—i.e. it should remain visually similar to 
its 2D equivalent even if mapped onto a complex curved 
surface. 


5.1.1. Mapping a 2D noise 


2D noise can be mapped onto surfaces through planar 
parametrization, exactly like regular texture maps. However, 


== 


Anisotropic noise 


` 


ensation 


Figure 8: Anisotropic filtering and compensation for para- 
metric distortions withanisotropic noise. (a) Anisotropic 
noise leads to higher image quality compared to isotropic 
filtering, as shown by the difference between close-ups. (b) 
Anisotropic noise compensates for parametric distortions 
to enforce a uniform noise aspect along the surface, as 
shown by the difference between close-ups. (Figure from 
[GZD08],©ACM, 2008.) 


this can introduce distortions and seams, breaking important 
properties of the noise such as uniform frequency content, 
continuity and whether the noise is band-pass. 


Goldberg et al.[GZD08] compensate for mapping distor- 
tions by locally adapting the noise content [Figure 8(b)]. This 
is a form of dynamic spectral control (Section 4.1), where 
the weights of the noise bands are driven to compensate for 
the distortions. The local distortion as well as its impact on 
the noise spectrum is estimated at every pixel. A noise with 
inversely pre-distorted frequency content is generated to ap- 
pear uniform along the surface. This is done by updating, in 
every pixel, the weights of the summed noise bands to ap- 
proximate the pre-distorted spectrum. This can be performed 
efficiently from a shader running on the GPU. This approach, 
however, only recovers from distortions and cannot hide the 
seams. Note that this idea was also suggested in the work of 
van Wijk [vW91, Figures 11 and 12]. 


5.1.2. Sampling a solid noise 


Noise can be applied onto surfaces by sampling a 3D noise 
function at every surface point. All noise functions are easily 
generalized to 3D and higher dimensions. Explicit noises, 
however, quickly induce a large memory cost because they 
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rely on pre-computed tables. The idea of sampling a 3D noise 
on surfaces was introduced by Perlin [Per85] and Peachy 
[Pea85]. This is often referred to as solid texturing. The 
approach, which popularized procedural textures, has several 
advantages: It is simple, memory consumption remains low, 
and the object appears as if carved out of a block of matter, an 
effect difficult to achieve otherwise. For a complete overview 
of solid texturing please refer to Dischler and Ghazanfarpour 
[DG01]. 


Cook and DeRose [CD05] observed that sampling a solid 
band-pass noise along a surface does not result in a band- 
pass noise on the surface. This is a consequence of the slice- 
projection theorem [Bra99, Mal93], which states that slicing 
in one domain corresponds to projection or integration in the 
other domain. Evaluating a 3D noise along a surface corre- 
sponds to slicing. Therefore, the power spectrum of the noise 
on the surface is given by integrating the band-pass power 
spectrum of the solid noise. However, this power spectrum 
is not band-pass anymore. Cook and DeRose additionally 
observed that the slice-projection theorem also provides a 
solution to this problem. Integrating a solid noise perpen- 
dicular to the surface corresponds to projection. Therefore, 
the power spectrum of the noise on the surface is given by 
slicing the band-pass power spectrum of the solid noise. This 
power spectrum is still band-pass. This provides a general 
method for obtaining a band-pass noise on a surface from a 
band-pass solid noise. 


5.1.3. Defining noise directly on the surface 


A last alternative is to define a noise directly on a surface, 
so that its features flow along the curvatures and naturally 
adapt to topology changes. This is difficult in general, but 
sparse convolution noises enable this approach: By locally 
splatting kernels the noise appears along the surface without 
having to resort to a global planar mapping. These ideas were 
hinted in earlier work [Cha07, 5.2] and further developed by 
Lagae et al. [LLDD09a]. In this latter work, the noise pattern 
is procedurally generated along a surface without any pre- 
processing such as computing a surface parameterization. At 
any evaluation point, only the 3D point coordinates and the 
surface normal are necessary to evaluate the 2D noise. For the 
case of anisotropic (oriented) textures, a direction field must 
also be provided to indicate the orientation of the texture. 
Several methods are available for the design of such fields 
(see, e.g. Fischer et al. [FSDHO7]). 


Solid texturing and surface noise produce different vi- 
sual effects. The first creates the unique feeling that the ob- 
ject is sculpted out of solid matter, whereas the second lets 
anisotropic textures ‘flow’ around the object. This is impor- 
tant when texturing, for instance, objects made out of fibres 
(straw basket, woven cloth, etc.). Figure 9 illustrates this 
idea. 


Figure 9: Difference in aspect of solid noise and surface 
noise. Straw hat textured with both a solid noise (left) and 
a surface noise (right). Left: The straw orientation is fixed 
in space, resulting in stretch on the side of the hat. Right: 
The straw orientation flows around the surface, producing 
the appropriate effect. There are no texture coordinates in 
both cases. (Figure from [LLDD09a], €ACM, 2009.) 


52. Filtering noise on surfaces 


An important consideration when mapping noise on surfaces 
is filtering of the frequency content when objects are seen 
at an angle or from a distance. This is crucial for rendering 
quality. Super-sampling is generally only necessary at geo- 
metric edges because textures are filtered, for instance using 
MIP-mapping. A major drawback of procedural textures is 
that such filtered lookups may not be available, requiring the 
use of super-sampling on the entire image. Because textures 
contain very fine details, and are seen from very close to far 
away, super-sampling will often not be able to solve the prob- 
lem entirely at reasonable cost. It is thus crucial to provide 
filtered sampling of procedural textures. 


In Appendix, we provide the necessary background to un- 
derstand filtering of signals mapped to surfaces. Although 
filtering is typically seen as a convolution in the spatial do- 
main, it can also be interpreted as a multiplication in the 
frequency domain. More specifically, the spectrum of the fil- 
tered noise is given by the multiplication of the spectrum of 
the unfiltered noise and the spectrum of the filter in texture 
space. The filter in texture space varies in each screen pixel 
because it is view-dependent. Figure 10 illustrates these 
concepts. 


We first describe how noise can be filtered (Section 5.2.1), 
and then discuss filtering of texture patterns obtained by 
applying transformations to noise values (Section 5.2.2). 


52.1. Filtering noise 


The key idea of filtering noise is to exploit the spectral con- 
trol offered by the noise to directly generate noise with the 
filtered power spectrum, rather than explicitly filtering unfil- 
tered noise. 


When noise patterns are obtained as a weighted sum of 
band-pass noises, a first approach is to cancel the contribution 
of bands whose frequency is too high. This approach is often 
referred to as frequency clamping [NRS82]. This works best 
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Figure 10: Anisotropic filtering with Gabor noise. Top: Un- 
filtered noise mapped on a tilted plane. The noise pattern is 
incorrect in the distance due to aliasing. Middle, from left to 
right: The power spectrum of the unfiltered noise, the filter 
for pixels in the red circle area, the power spectrum of the 
filtered noise for these pixels. This last spectrum is simply 
the product, in the frequency domain, of the filter and the 
noise power spectrum. Bottom: Same noise but properly fil- 
tered. Aliasing is entirely removed. (Figure from [LLDD09a], 
©ACM, 2009.) 


if the noise is narrowly band-pass (i.e. the ring in the spectrum 
has to be thin and well defined). Perlin noise [Per85] is only 
weakly band-pass, making frequency clamping difficult to 
tune. Cook and DeRose alleviate this issue by providing a 
noise with better defined band-limits [CD05]. 


Both of these noises, however, are isotropic and the clamp- 
ing cannot account for the anisotropy of the filter. On tilted 
surfaces one must compromise between over-blurring or 
residual aliasing. Goldberg et al. [GZD08] obtain higher 
quality filtering since their pre-computed noise bands are 
oriented. Each band corresponds to a noise pattern with lim- 
ited frequency content along a given orientation (Figure 11). 
By adapting the weights of the oriented bands with respect to 
the anisotropic filter, the noise content adapts to non-uniform 
perspective distortions [Figure 8(a)]. 


Lagae et al. [LLDD09a] exploit a unique property of their 
noise. The noise is obtained as a sum of Gabor kernels. Each 
Gabor kernel corresponds to a Gaussian in the frequency 
domain. It is possible to filter each individual kernel by com- 
puting the product, in the frequency domain, between the 
Gabor Gaussian and the filter Gaussian. Because Gaussians 
are closed under multiplication, the product is a third Gaus- 
sian. This new Gaussian can be interpreted as a filtered Gabor 
kernel. The parameters of this new kernel are used instead 
of those of the original, unfiltered kernel. This directly gen- 
erates a noise with a filtered spectrum. Contrary to previous 
methods exploiting a discretization of the spectrum in distinct 
bands, this approach allows analytical filtering. 


1. Uniform white noise 


2. Frequency domain decomposition 


3. Inverse transform 


Figure 11: Anisotropic noise generation. Illustration of 
spectral noise generation. The frequency domain decompo- 
sition has three orientations. Three oriented subbands at 
the same scale and their corresponding spatial domain im- 
ages are shown, which are stored as textures. (Figure from 
[GZD08], ©ACM, 2008.) 


5.2.2. Filtering noise-based procedural textures 


Noise patterns are rarely used directly to produce textures. 
Patterns are generated by applying several functions to the 
noise, such as absolute values or sine waves. In addition, the 
noise is often coloured by remapping its values to a piecewise 
linear colour ramp [EMP*02]. Figure 12 illustrates how a 
marble texture is built from a solid noise. 


Because most of these additional operations are non-linear, 
starting from a filtered noise value does not guarantee that the 
resulting texture is also filtered. Although this approximation 
is acceptable when the function applied to the noise is very 
smooth, proper filtering is in general necessary. For example, 
consider a black and white pattern obtained by thresholding 
the noise. A correctly filtered version should progressively 
introduce blur in the transition areas. However, when only 
filtering the noise, the transitions will remain sharp due to 
the subsequent thresholding. 


Although this problem remains unsolved in the general 
case, several approaches provide good approximations when 
the operations applied to the noise can be summarized in a 
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Figure 12: Procedural texture creation. The marble vase 
(left) is obtained from two components: A colour map re- 
peated along the x-direction in space (middle), perturbed by 
a solid noise (right). The final colour is obtained as C(x 
+ N(x, y, z)) where C is the 1D colour map, N the noise 
and x, y, z the surface point coordinates. (Figure based on 
[LLDD09a],OACM, 2009.) 


1D colour table. The final colour is obtained as C(N(u)), 
where N(u) is the noise value at u and C the colour table. 


Rhoades et al. [RTB*92] filter the colour table C rather 
than the noise function. They return the average colour 
over a small interval [N(u) — ô, N(u) + 4]. This is con- 
veniently evaluated using MIP-mapping on the 1D colour 
table. The size of the interval ó is computed from the filter 
size and the maximum gradient of C with respect to u. Hart 
et al.[HCK99] further refine this approach using the local 
noise gradient (most noises are differentiable, either through 
finite differencing or analytically). Although this works well 
in many cases, one source of error is that the noise value 
and gradient are evaluated on the unfiltered noise. Lagae 
et al.[LLDD09a] rely on a similar mechanism. They achieve 
accurate filtering by estimating the noise value range ó from 
the loss in noise variance due to filtering. This is only possi- 
ble because an analytical expression of the noise variance is 
available. 


Several authors have investigated more general methods 
for filtering procedural textures. Heidrich et al. [HSS98] pre- 
sented a method to obtain an average value of a procedural 
shader with an error bound over a finite area using affine 
arithmetic. Olano et al. [OKS03] presented a method for au- 
tomatic shader level-of-detail using an automatic system for 
shader simplification. 


6. Analysis of Procedural Noise Functions 


In this section, we give a detailed analysis of procedural 
noise functions. We introduce analysis tools (Section 6.1) 
and present analysis results (Section 6.2). 


6.1. Analysis tools 


Motivated by our definition of noise as a stationary and 
normal stochastic process, we introduce analysis tools for 


estimating the power spectrum and the amplitude distribu- 
tion of a noise function. The estimated statistics of a noise 
function can provide insight into the noise function, also in 
the case when the expected statistics are known. For example, 
differences in expected and estimated statistics might reveal 
mistakes in implementation. 


6.1.1. Power spectrum 


We estimate the power spectrum of a noise function using 
Bartlett’s method of averaging periodograms [Bar78]. The 
periodogram is a simple estimator for the power spectrum, 
defined as the magnitude squared of the Fourier transform 
[PVTF02, 13.4]. However, the periodogram is very noisy. 
This is because the periodogram is a white noise process 
with as mean the power spectrum [PP02, 12.2]. Averaging 
periodograms of different instances of noise results in a less 
noisy estimate for the power spectrum of a noise function. 
We inspect both the power spectrum estimate as well as the 
periodogram, because averaging periodograms averages out 
noise but can also average out features. We radially aver- 
age the power spectrum of an isotropic noise function, be- 
cause the power spectrum of an isotropic noise function is 
radially symmetric. Note that these methods are also used 
for power spectrum estimation of Poisson disk distributions 
[U1i88, LDOS]. 


6.1.2. Amplitude distribution 


We estimate the amplitude distribution of a noise function us- 
ing a histogram of noise values. We plot a Gaussian function 
with the expected variance, in case this is known for the noise 
function, or with the estimated variance, as a reference. 


6.1.3. Other analysis tools 


Yoon et al. [YLC04, YLO08] presented several other analysis 
tools for measuring the quality of a noise within an optimiza- 
tion procedure. Yoon et al. used a Chi-square goodness-of-fit 
test to measure the quality of the amplitude distribution of 
the noise, a test based on autocorrelation to detect periodicity 
in the noise, and a band-pass test to measure how band-pass 
the noise is. The band-pass test is inspired by wavelet noise 
[CD05], and is based on the difference between the noise 
and a down-up-sampled version of the noise. However, for 
our purpose, these analysis tools are subsumed by the ones 
above. 


6.2. Analysis results 


We have analysed Perlin noise (using the implementation 
of [Per02], see Section 3.1.1), sparse convolution noise (as 
presented in [Lew89], see Section 3.3.1), wavelet noise (as 
presented in [CD05], see Section 3.2.1), anisotropic noise 
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(as presented in [GZD08], see Section 3.2.2), better gradient 
noise (as presented in [KKS08], see Section 3.1.2) and Gabor 
noise (as presented in [LLDD09a], see Section 3.3.3). 


The parameters of the noise functions were selected to pro- 
duce a noise with a principal frequency of 1/32 for the spatial 
domain and 1/4 for the frequency domain. Note that both do- 
mains require different parameters for optimal visibility. The 
spatial domain images were tone-mapped by linearly map- 
ping a range of three standard deviations to intensity, and the 
frequency domain images by linearly mapping the expected 
maximum of the power spectrum to an intensity value of 
80%. One hundred periodograms were used to compute the 
power spectrum estimate. 


We present the results of our analysis in Figure 7. 


In row 1, we show the noise generated by the noise func- 
tions. Wavelet noise, anisotropic noise, better gradient noise 
and Gabor noise have a very similar aspect. Sparse convolu- 
tion noise has a different aspect, because it is not designed 
to be band-pass. Perlin noise has a slightly different aspect, 
because the noise is zero at every integer lattice point, and be- 
cause of an undesired axis-aligned anisotropy. Wavelet noise 
has a very subtle different aspect, because of an undesired 
axis-aligned anisotropy. 


In row 2, we show the amplitude distribution of the noise 
functions. All noise functions except Perlin noise have an 
approximately Gaussian amplitude distribution. The ampli- 
tude distribution of Perlin noise contains undesired arte- 
facts, because of the limited number of random gradient 
vectors, and because the noise is zero at every integer lattice 
point. 


In rows 3-5, we show the periodogram, the power spec- 
trum estimate, and the radially averaged power spectrum of 
the noise functions. All noise functions except Perlin noise 
and sparse convolution noise are approximately band-pass. 
Perlin noise is only weakly band-pass, which might lead to 
problems with aliasing and detail loss [CD05]. Sparse convo- 
lution noise is not designed to be band-pass. The horizontal 
features in the periodogram of Perlin noise are caused by 
an undesired correlation in the hash function [KKS08]. The 
axis-aligned square feature in the power spectrum of wavelet 
noise is caused by the separable B-spline [CD05]. Both fea- 
tures indicate an undesired axis-aligned anisotropy in the 
noise. Note that the horizontal features in the periodogram 
of Perlin noise are much less visible in the power spectrum, 
which is an example of a case where averaging periodograms 
can also average out features. 


We conclude from our analysis that the noise functions are 
often very different in terms of visual aspect, power spectrum 
and amplitude distribution. However, as we will show in 
Section 7, every noise function represents a specific trade-off 
between a set of features, and this analysis only takes into 
account a small part of this set of features. 


7. Comparison of Procedural Noise Functions 


In this section, we give a detailed comparison of procedural 
noise functions, based on the previous sections. 


We compare the same noise functions as the ones we have 
analysed in Section 6. We present the results of our com- 
parison in Table 1. It is important to note that several de- 
velopments presented in later methods are also applicable to 
earlier methods. For example, several developments in better 
gradient noise and Gabor noise are applicable to Perlin noise 
and sparse convolution noise respectively. In the table, we 
compare the methods as presented in the cited works, while 
in the discussion, we generalize. 


In part (a) of the table, we compare to which degree the 
noise functions adhere to the definition of procedural noise 
function (Section 2.3). Storage requirements and periodic- 
ity are generally linked. Explicit noises have high storage 
requirements, whereas other noises have low storage re- 
quirements. Several authors presented methods to improve 
storage requirements and periodicity, for example, noise 
tiles [CD05, YL08, Lag09] and long-period hash functions 
[LD06]. Sparse convolution noises are non-periodic and have 
minimal storage requirements. 


In part (b) of the table, we compare the noise functions in 
terms of modelling (Section 4). Band-pass noise functions 
achieve spectral control using a weighted sum of noise oc- 
taves, while sparse convolution noises achieve spectral con- 
trol using the kernel. 


In part (c) of the table, we compare the noise functions in 
terms of noise on surfaces (Section 5). All noise functions 
generalize to arbitrary dimensions, and can therefore support 
solid noise, although the high storage requirements of explicit 
noises can be problematic for solid noise. Sparse convolution 
noises can support surface noise, and all noise functions that 
support band-pass solid noise can support band-pass surface 
noise. 


In part (d) of the table, we summarize the noise func- 
tions in terms of filtering (Section 5). All noise functions 
that are band-pass can support isotropic filtering, while only 
noise functions that support anisotropic noise can support 
anisotropic filtering. 


In part (e) of the Table, we summarize the analysis of 
the noise functions (Section 6). Perlin noise does not have 
a Gaussian amplitude distribution and is only weakly band- 
pass. Sparse convolution noises are band-pass when the ker- 
nel is band-pass. 


In part (f) of the Table, we compare the noise func- 
tions in terms of speed. In our experience, sparse con- 
volution noises, wavelet noise and better gradient noise 
are generally slower than Perlin noise, whereas anisotropic 
noise is generally faster. Note that this is a subjective com- 
parison. Sparse convolution noises offer a speed versus 
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Table 1: Comparison of procedural noise functions. 


Perlin Sparse convo- Wavelet Anisotropic Better gra- Gabor 
noise lution noise noise noise dient noise noise 
[Per02] [Lew89] [CD05] [GZD08] [KKS08] [LLDD09 a] 
Category Lattice gradient Sparse convolution Explicit Explicit Lattice gradient sparse convo. 
(a) Definition—procedural 
Storage requirements! O(N) O(1) O(N?) O(N?) O(dN) O(1) 
Contin, no discrete data NA K af: n4 af 
Non-periodic? J J 
Parameters? Weights Kernel Weights Weights (aniso.) Weights Kernel param. 


(b) Modelling—spectral control 
Spectral control 
Anisotropic noise 

(c) Surfaces—noise on surfaces 
Solid noise x af: 
Surface noise 
Setup-free surface noise 

(d) Surfaces—filtering 
Isotropic filtering 
Anisotropic filtering 

(e) Analysis 
Band-pass power spectrum 
Gaussian amplitude distribution WA 

(f) Comparison—speed 
Speed* 0 — 
Quality/speed tradeoff NA 


Weight. sum Kernel 


Weight. sum Weight. sum (aniso. Weight. sum Kernel param. 


" J 
" 4 J 
y ", Y Y 
H d Y 
v Y y 
J v 

y y " J 
y y v Y 
T + = pes 
J 


Note: Comparison of Perlin noise, sparse convolution noise, wavelet noise, anisotropic noise, better gradient noise and Gabor noise. See 
Section 7 for a correct interpretation of this comparison. (Table based on [LLDD09a], €ACM, 2009.) 

! Storage requirements are expressed in function of the period N and the number of dimensions d. 

? Non-periodicity does not take into account the inherent periodicity of computer calculations. 


? Parameters that are required for spectral control. 
^ Speed is expressed relative to the speed of Perlin noise. 


quality tradeoff, but remain slower for acceptable quality 
levels. 


We conclude from our comparison that every noise func- 
tion represents a specific trade-off between a set of features, 
and that the noise function that is best suited for a specific 
application depends heavily on the requirements of that par- 
ticular application. 


8. Conclusion 


In this survey, we have provided a critical view on proce- 
dural noise methods, including recent solutions developed 
in the last 8 years. We started by providing a well-founded 
definition of noise, with a theoretical grounding in stochastic 
processes. These definitions allowed us to provide a unified 
classification of the most important procedural noise solu- 
tions, providing the reader with a coherent view of the field. 
In particular, we have classified procedural noise solutions 
into lattice gradient, explicit and sparse convolution noises. 
We underline the importance of spectral control when mod- 


elling patterns and visual detail with noise. In many cases, 
noise is applied to surfaces: we have distinguished the dif- 
ferent ways to do this, notably via mapping of 2D noise, 
solid noise and local kernel splatting. Several important is- 
sues arise when applying noise to surfaces, most of which 
have been treated in recent work. In particular, these relate 
to filtering of noise, as well as the procedural textures based 
on noise. 


We have used the power spectrum and amplitude distribu- 
tion to analyse procedural noises. This analysis helps explain 
some of the difference between noises and, in some cases the 
differences in the resulting visual aspect. We also provided 
a comparison of the various noises, based on the features 
that each solution provides. These include storage require- 
ments, the way the power spectrum is controlled, whether 
anisotropic noise is provided, how each noise can be applied 
to surfaces and consequent filtering solutions, and of course 
speed of computation. The conclusion of this comparison is 
that each solution presents a different tradeoff. Each appli- 
cation needs to determine the relevant importance of each 
feature to determine which noise is most appropriate for the 
problem at hand. 
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8.1. Future work 


We next discuss several challenging directions for future 
work. 


8.1.1. Fine-grained control over the power spectrum 


In Section 2, we have explained that noise is completely char- 
acterized by its power spectrum, and in Section 4, we have 
discussed several methods for controlling the power spectrum 
of noise. However, these methods only offer coarse-grained 
control over the power spectrum. Fine-grained control over 
the power spectrum would bridge the gap between noise and 
stochastic texture. This is illustrated by Lewis [Lew86, Figure 
6] and by Lagae et al. [LLDD09a, Figure 7]. Finer control 
over the power spectrum is an interesting direction for fu- 
ture work; it is particularly important for the development of 
noise-by-example methods. 


8.1.2. Non-Gaussian and non-stationary stochastic 
processes 


In Section 2, we have defined noise as a stationary and 
Gaussian stochastic process. Generalizing this definition di- 
rectly suggests two ways to extend noise: non-stationary and 
non-Gaussian stochastic processes. Non-stationary stochas- 
tic processes correspond to spatially varying noise. Lewis 
[Lew86] already hinted at spatially varying sparse convo- 
lution noise, van Wijk [vW91] used spatially varying spot 
noise for visualizing scalar fields, and Lagae et al. [LLD10] 
recently explored spatially varying Gabor noise. However, all 
the above solutions have limitations. Developing a general 
and efficient solution allowing spatial variation of all noise 
parameters is an interesting direction for future research. 
Non-Gaussian stochastic processes correspond to more gen- 
eral random patterns. This could bridge the gap between 
noise and non-stochastic texture, or between noise and para- 
metric texture synthesis [PNNT96, PSOO]. Although some 
authors have investigated non-Gaussian processes, for ex- 
ample, Lewis [Lew86], Gagalowicz and Ma [GM85], and, 
more recently, Chainais [Cha07], research into non-Gaussian 
noise is very limited. The development of non-Gaussian noise 
would provide a powerful tool to model much richer patterns. 
It is however unclear whether such an approach is the most 
effective way to obtain such results. In addition, developing 
general solutions for non-Gaussian noise requires substantial 
future research. 


8.1.3. Understanding and controlling phase 


In Section 2, we also noted that a noise is fully specified 
by its power spectrum, but has random phase. This charac- 
terization suggests that the distinction between a stochastic 
texture and a structured texture or pattern is manifested in the 
non-random phase of the latter. Identifying and controlling 


useful information in the phase might be another approach to 
bridging the gap between stochastic and structured textures. 
This is an unexplored topic. 


8.1.4. Faster noise and benchmarking noise 


Because the introduction of Perlin noise more than two 
decades ago, making noise faster continues to be impor- 
tant, especially for industrial applications. This is perhaps 
best illustrated by the famous quote: 90% of 3D rendering 
time is spent in shading, and much of that time is spent com- 
puting Perlin noise). In Section 7, we have not performed 
an extensive benchmark to compare the speed of different 
noise functions. This is mainly because of the practical is- 
sues involved, such as coding optimized implementations for 
different architectures (CPU, GPU). However, such a bench- 
mark would be interesting. As mentioned in the introduc- 
tion, the relative cost of computation and memory access has 
changed significantly over time, to the point where current 
program execution times are often dominated by memory 
access. This is one reason for the renewed interest in pro- 
cedural noise, that is the hope of exchanging bandwidth for 
computation. 


8.1.5. Better authoring tools 


In Section 4, we have mentioned some authoring tools for 
noise. However, the work done in graphical user interfaces 
for noise, in noise editing, and in noise by example is very 
limited. 


8.1.6. Filtering procedural textures based on noise 


In Section 5, we have discussed filtering of procedural tex- 
tures based on noise. Although solutions for specific cases are 
available, the general problem of filtering procedural textures 
based on noise is still unsolved. 
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Appendix. Background on Texture Filtering 


This section explains the basic concepts for filtering signals 
mapped to surfaces. We describe here the general approach, 
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which applies both to bitmap textures, procedural textures or 
any other signal applied to a surface. The intent is only to 
introduce the basics concepts needed for Section 5.2. Please 
refer to Heckbert [Hec89] for a more detailed account of 
texture filtering. 


The filter to be applied to the texture varies in each screen 
pixel, as it depends on the viewing condition. In screen space, 
each pixel is typically modelled as a small Gaussian filter 
(but this could be any other filter). The colour of the pixel 
should be obtained as the integral of the filter multiplied 
by the visible portion of the texture. In other words, this 
evaluates at the pixel the convolution between the filter and 
the projected texture. To simplify computations, the filter 
is back-projected in the texture domain through the view 
transformation and through the mapping used to apply the 
texture onto the visible surface. This results in a distorted 
Gaussian in the texture domain (Figure 10). 


More precisely, the filter for pixel (x, y) in image space is 
the Gaussian: 


fay= ga (A.1) 


1 
OoN 2T 
where o is the width of the Gaussian in pixels (typically 
o = 1.0). The corresponding filter in the texture space is the 
Gaussian 


1 ene vo ye v]T 
o 


o2z ) 
where J is a local affine approximation (Jacobian) of the 
mapping from image to texture coordinates. It is easily esti- 
mated in every pixel using the screen space derivatives of the 
texture coordinates (rasterization) or ray differentials (ray- 
tracing): 


fC uva 


(A.2) 


du du 
dx d 
J= i ; (A.3) 
dv dv 
dx dy 


Note that this Gaussian filter remains a Gaussian in the spec- 
tral domain. The corresponding frequency domain filter is 
the Gaussian: 


—22?c? JT T 
FE(II ful") = e207 fu fold ITU fol" (A 4) 


We thus can easily relate the filter to the spectrum of the 
texture. A properly filtered texture will have for spectrum the 
multiplication between its unfiltered spectrum and the filter 
spectrum. This follows from the convolution theorem: The 
spatial domain convolution simply becomes a multiplication 
in the spectral domain. 
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